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Nonlinear coupling between scissors modes of a Bose-Einstein condensate 
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We explore the nonlinear coupling of the three scissors modes of an anisotropic Bose-Einstein con- 
densate. We show that only when the frequency of one of the scissors modes is twice the frequency of 
another scissors mode, these two modes can be resonantly coupled and a down conversion can occur. 
We perform the calculation variationally using a gaussian trial wave function. This enables us to obtain 
simple analytical results that describe the oscillation and resonance behaviour of the two coupled modes. 
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I. INTRODUCTION 

Similar to monopole and quadrupole breathing modes of a gaseous Bose-Einstein condensate [1-7], scissors modes 
were first studied theoretically [8,9] and subsequently observed experimentally [10]. The scissors modes are, however. 
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lO i rather special since they directly manifest the superfluid behavior of these atomic gases. Moreover, the recent exper- 
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imental studies appear to show a resonance behavior between two coupled scissors modes [11]. From a theoretical 
point of view this is interesting because a linear-response calculation can neither account for the coupling nor for 
the resonance behavior [8,9]. Therefore, a first step towards an explanation of these experimental observations is 
to perform a calculation that goes beyond linear-response theory and accurately takes into account the mean-field 
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^ interaction that couples the scissors and breathing modes. In this paper we present a simple variational method for 



calculating the frequencies of these various modes and their couplings beyond the linear-response. We perform our 
calculation at zero temperature and therefore do not consider the damping rates of the scissors modes [12,13]. 

The main idea behind our method is to use a time-dependent gaussian ansatz for the groundstate wave function to 
derive the equations of motion of the breathing modes and the scissors modes. Then we expand the resulting equations 
of motion in deviations from equilibrium. In first order, i.e., linear-response, we recover the expected uncoupled set of 
equations [8,9]. The second-order calculation produces a set of coupled equations which show that we need to excite 
all three scissors modes in order to get a nonzero coupling. At higher orders we, however, find that we only need 
to excite two modes to get coupling. Furthermore, we actually find under certain conditions a resonance behavior 
between these two modes. 

The layout of the paper is as follows. First, we rederive in Sec. II the frequencies of the scissors modes in the 
linear-response limit. In Sec. Ill we extend the calculation first to second, and then also to higher orders, which 
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ultimately lead to a resonant coupling. In Sec. IV we solve the equations of motion analytically near the resonance 
using an envelope function approach. In Sec. V we end with a discussion of our results. 

II. FREQUENCIES OF THE SCISSORS MODES 

We start by considering a Bose-Einstein condensate trapped by the following harmonic potential 



V{r)^-m{ulx^+u,y+ulz^), (1) 



where w^,, ujy, and Wz are the angular frequencies of the trap, and m is the atomic mass. A scissors mode in a 
Bose-Einstein condensate is associated with an irrotational flow with a velocity field of the form v(r) ex V(xy), if 
the motion is taking place in the xy-plane [8]. Similar expressions hold for the two other Cartisian planes. These 
kind of modes can be excited by a sudden rotation of the equilibrium axes of the trap. To such a perturbation the 
condensate will respond by oscillating around the new equilibrium axes. For example, to excite a scissors mode in 
the xy-plane, we rotate the x and y-axes of the trap slightly around the z-axis. If the angle of rotation is sufficiently 
small, the scissors mode can be approximated by a simple oscillation of the condensate around the new equilibrium 
axes. On the other hand, if the axes change through a large angle this method excites the to = 2 quadrupole mode, 
where to. labels the projection of the angular momentum along the axis of symmetry. The maximum angle for which 
the scissors mode is defined increases with deformation of the trap [14]. 

To account for all three scissors modes in the three Cartisian planes we employ the following trial function for the 
condensate order parameter 

■i/'(r, t) = A(t) exp (-fea;(i)x^ - by{t)y^ - bz(t)z^ - Cxy{t)xy - Cxz{t)xz ~ Cyz{t)yz) , (2) 

where hi and Cij , are complex time-dependent variational parameters and 

y/N2^/^ J ^ 

— 3/4 \/ ^xy,r^xz.rCyZjr ~r ^Ox^r'^y^r'^z^r \^z^r^xy,r ' ^y^^^xz^r "*" ^x,rCyz^r ) ' V*^/ 

This value of the prefactor A{t) guarantees the normalization of the square of the wave function ip{r, t) to the 
total number of condensed atoms N. Here hi^r and Cy,, are the real parts of hi and c^, respectively. The first 
set of parameters, hi, give rise to the well-studied breathing modes which, for axially symmetric traps, are called 
the monopole and quadrupole modes depending on the value of m being equal to zero or two, respectively. The 



parameters c^ on the other hand determine the three scissors modes. The equations of motion for these variational 
parameters can be derived from the lagrangian 

where E[^,ip*] is the usual Gross-Pitaevskii energy functional given by 
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E[ip,^p*] = / dr 



Zm 2 



(5) 



Here T^^ is the two-body T-matrix, which for the atomic Bosc-Einstcin condensates of interest is related to the 
s-wave scattering length a through T^^ = Airafi^ /m. 

Inserting our trial wave function into the lagrangian and scaling frequencies with Hi ~ {iOxLOyLOz)^ and lengths with 



a = yjh/rau}, it takes the dimcnsionlcss form 

L[6, c]/N = {oixbxA + ctyby,i + '=^zbz,i] IQ (6) 



- - {0lx{^\hx'^ + \Cxy'^ + \Cxzf) 
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axi^l + OLyU?y + OL^u?^ IQ 



where Q — 271^ A^ /N"^ , ax — 46j, ,.&z,r — Cy^.r, cty — 4:bx.rbz,r — clcz,ri '^z — 4:bx.rby^r — c^j/,ri a'^d the dot corresponds 
to a time derivative. In addition 7 — Na/a is the dimensionless parameter that represents the strength of the mean- 
field interaction. Minimizing the lagrangian with respect to the 12 variational parameters, we get a set of 12 coupled 
equations of motion. The resulting equations of motion are rather lengthy and complicated. A significant simplification 
takes place if we expand these equations in the deviation of the variational parameters from their equilibrium values, 
i.e., in Sbi{t) = cxp (—iujt)(bi — b\ ) — cxp(—iujt)5bi and Scij{t) = cxp {—iujt)(cij — c^.- ) = cxp{—iujt)Scij, where 
uj is the yet unknown eigenfrequency of the modes, and bl and q- denote the equilibrium values of bi and Cy, 
respectively. These latter values can simply be obtained by setting the time derivatives in the equations of motion to 
zero. For large condensates in the so-called Thomas-Fermi regime we can ignore contributions from the kinetic energy 
[15] and the equilibrium variational parameters take the simple form 



cl? = 0. (7) 

It is required that 7 3> 1 for the last equation to be valid. To first order in the deviations the equations of motion 
read simply 

M-P = 0, (8) 

where the vector P = {Sbx,r, Sby^r, Sb^^r, Sb^^i, 6by^i, 6bz,i, dc^y^r, Scxz,r, 6cyz,r, 5cxy,i, Scx^a, Scyz,i) contains all the pos- 
sible fluctuations, and the matrix M is given by 

/ Inbreathing q \ 

^^^ I A lyrscissors I i W/ 

where Mbroathmg g^^^ jyi-scissors g^j.g gjygjj^ explicitly in Appendix A. It is clear from the last equation that to linear 
order the breathing modes and scissors modes are uncoupled. The dispersion relation of these modes can be obtained 
by setting the determinant of M to zero. This results in 

X {uj^ - SloIlo^ + 8w^a)2 - 20uj^) ^ 0, (10) 



where il^y = ^aj^+u^, fl^^ = y/uj^+uj^, ny^ = ,JZf+Zj, uj^ ^ uj^ + uj^ + uj^ , uj^ ^ uj^ujI + ojIujI + w^w^, and 
Lo^ = [uJxUJyOJz)"^ ■ The zeros of the first three factors in the left-hand side of the last equation give the frequencies of 
the scissors modes, whereas the zeros of the sixth-order polynomial give the frequencies of the breathing modes. 

So far we have not produced a new result, since these frequencies have been calculated previously [8,9,16]. However, 
we have demonstrated a new way of calculating the frequencies of scissors modes in a relatively simple manner. The 
most important part of this paper is therefore contained in the next section, where we consider also the nonlinear 
effects produced by the mean-field interaction. 

III. BEYOND LINEAR RESPONSE 

In this section we consider the equations of motion for the variational parameters by taking into account several 
higher-order terms in the deviations 5bi and 5cij. We have calculated these equations analytically up to second order, 



but they turn out to be rather lengthy and contain terms that couple the breathing modes and the scissors modes. 
As we discuss below, in the present experiments with axially symmetric traps the coupling of the scissors modes with 
the quadrupole mode is always of importance, but we leave the treatment of this more complicated situation to future 
work. For simplicity, therefore, we focus here on triaxial traps, in which case we can ignore the breathing modes. 
This is achieved by simply setting Sbi = in the full equations of motion. 

Up to second order, the remaining 6 first-order equations for Scxy^r ■ ■ ■ Scyz,i can be reduced to 3 second-order 
equations for Sc^y^r, Sc^z.r, and Scyz.r by eliminating the imaginary parts. We find in detail 



'J'^xy,r ^ ^''xy^'^xy.r ^ \ ) 1 I 4 I 



2/5 



-nljCxz.r5Cyz,r + 2 f J^"^^ 1 5. 



3^ 2/ 



yz 



(^xz.r^^yz,r 



(11) 



Scxz,r + nljCxz,r + {-y/'l'^' ('^) 



2/5 



2;^^*-'3:y,r^Cyj2,7- -^r ^ 



and 



5Cy.,r + nljCyz,r^{^f/W"{^-^) 



2/5 



^ty^OCxy^r^^Cxz^r i ^ 



i Lxy^ ^yz 



i Lxy^ I'XZ 



^(^xy.r^^yz. 



y.r'^'~-yz.r 



(12) 



^Cxy,r^Cxz,r 



(13) 



These equations show that if we ignore the second-order terms we recover our previous three uncoupled scissors modes 
with frequencies flxy, ^xz and flyz. It is interesting to observe that the coupling terms couple the three scissors modes 
such that if one mode is not excited the other two modes will not be coupled. In other words, to get a nonzero coupling 
up to this order, the three modes must be excited simultaneously. We believe that this important result is not an 
artefact of the gaussian approximation, but also holds for an exact calculation using the Bogoliubov theory. 

For higher-order couplings the last conclusion is no longer true. Two modes can then be coupled, even when the 
third is not excited. From now on, therefore, we assume without loss of generality that only the Sc^y and Sc^z modes 
are excited, while Scyz — 0. Moreover, to investigate the possibility of a resonant coupling between these two modes, 
we have considered coupling terms up to ninth order. Similar to the second-order case explained above, the equations 
of motion can thus be expressed as 



^^xy.r ~r ^ ^xy^^xy,r i / ^ 

j+k+l+m<9 



^jklmy^Cxy.r ) y^^xy.r) y^^xz^r) [y^xz^rj 



(14) 



^C-xz^r \ i L^^OCxz,r ~r 



Z' Pjkhn \P(^xz.r) \^(^xz.r ) x^^xy.r } x^^xy.r } j 

j+/c+/+m<9 



(15) 



where j,k,l,m ^ 0, 1, 2, 3 and the sum j + k + l + m does not exceed 9, which is the order up to which we have chosen 
to expand the equations of motion. The coefficients ajUm, (ijkim sue given in terms of the trap parameters lu^, Wy, 
Wz, and the interaction parameter 7. A resonance between these two modes takes place when the frequency of the 
coupKng terms is equal to the frequency of the zeroth-order term, i.e., the first two terms. The frequency of each 
coupling term in the above summations is determined by substituting for 5cij and 6cij their zeroth-order solutions. 
Therefore, the latter frequencies will be a certain linear combination of the zeroth-order scissors modes frequencies. 
Imposing the above resonance condition on each coupling term thus results in a relation between the two scissors 
modes frequencies fl^y and Qxz- Inspecting all coupling terms up to the ninth order, we found a very small number 
of terms for which the relation between Q^y and Qxz can be satisfied by real values of cux, ojy, and ujz- Ultimately, we 
find only a resonance when either 

or 

is satisfied. The first resonance condition leads in the experimentally relevant axially symmetric case, where lUx = 
ujy ~ ujz/\/\, to a value for the anisotropy ratio A that is equal to 1. This clearly corresponds to a spherically 
symmetric condensate. Since in this case the scissors modes are degenerate with the quadrupolc breathing modes, 
which we have ignored here, we focus from now on only on the second resonance condition. 

It is interesting to mention that the second resonance condition is exactly the same as the one observed experimen- 
tally [11]. For a resonance of this kind the resonant coupling terms turn out to be of seventh order, and the equations 
of motion in that case read 

dCxy,r + i^ljCxy^r + P {5Cxy^rT {5Cxy^rf [SCxz.rT + V {SCxy^r)^{SCxz,r)'^ = 0, (18) 

6Cxz,r + ^lz^Cxz,r + a {SCxz.rf {SCxy-rf {SCxy^r)"^ = 0, (19) 

if we neglect all nonresonant terms. Here a, (3 and "q are functions of w^;, tOy, lOz, and 7 that are given explicitly in 
Appendix A. It is important to note here that our neglect of the nonresonant terms is justified when we are close 
enough to resonance. This is similar to the rotating-wave approximation known from quantum optics. We see that 



the coupling terms indeed lead to the above-mentioned resonance condition, by inserting in them the zcroth-order 
solutions, i.e., Scxy^r oc cicp {—iQ^yt) and 5cxz,r cc cxp (iQ^zt)- For example, the coupling terms in Eq. (18) have 
a total frequency of 3^xz — -^^xy Separating out Scxy.r oc exp (—iflxzt) as a prefactor for the whole equation, the 
coupling term oscillates thus as 2 A = Aflxy — 2ilxz- Therefore, a resonance takes place when A = 0, i.e., when the 
condition in Eq. (17) is met. Similarly, the coupling term of Eq. (19) is also oscillating with a frequency of 2A. 

IV. SOLUTION OF THE EQUATIONS OF MOTION NEAR RESONANCE 

Sufficiently close to resonance we can write Scxy and Scxz as a product of two functions. One of them describes the 
slow envelope and the other the fast oscillation with the uncoupled scissors-mode frequency. In particular, we have 

Scxyit) = g(t) exp {iQ,xyt), (20) 

and 

Scxz{t) = , fit) exp {-mxzt), (21) 

where g(t) and f{t) are the slowly varying envelope functions. Substituting these expressions into Eqs. (18) and (19), 
ignoring second-order time derivatives of f{t) and g{t), and then eliminating g(t), we obtain the following equation 
for fit) 

-(^/ + 2A/)/ + ^£/2 = 0, (22) 

where 

" ~ ^^y • ^ ^ 

This equation has a solution of the form 

fit) ^ [Ci + C2 exp (2i Ai)] T^ , (24) 

where Ci and C2 are two constants of integration that are determined by the initial conditions. Note that the relevant 
quantity here is |/(i)|, which represents the actual envelope of the oscillation and is given by 

\fit)\ = [d + Cl + 2C1C2 cos (2Ai)] ^^ . (25) 



In first instance we might think that the real part of /(t) is the relevant quantity. However, in Eqs. (20) and (21) 
we should in principle have taken the real part of the right-hand side. If we do that we automatically are lead to the 
condition that |/(i)| is the envelope of the oscillation. 

For definiteness sake let us take the initial conditions /(O) = fr{0) = ./o and /(O) = /i(0) = /o, where fr{t) and 
fi{t) are the real and imaginary parts of /(i), respectively. Physically, this set of initial conditions corresponds to 
exciting the scissors modes in the xz and the xy-planes simultaneously. This should be performed experimentally 
by initially rotating the condensate in the a;z-plane by an angle 60 and around the 2;-axis by an angle (J)q and then 
releasing the condensate. The initial angles 9q and 00 are related to the constants /o and /o by 

/o = 2|6f -5(")|cos0o sin^o, (26) 

/o = ^(2|6(")-6(")|cos0o sin0o)->i")-&f Icos^o sin^o- (27) 

With these initial conditions the constants Ci and C2 are given by 

gi = .r- + ^°-^%^ ' '\ (28) 

g2 = - -^°-^%^~"^ (29) 

Using these expressions and the experimental parameters from Ref. [10], we give in Fig. 1 the real part of Scxz{t). 
This clearly shows how the energy is being exchanged between the two modes. 

An interesting property of Eq. (22) is that exactly on resonance, i.e., A = 0, its solution becomes nonoscillatory. 
Indeed we find that in this case the solution is 

f{t) = {Cs + C^t)T^, (30) 

where again C3 and C4 are constants that are determined by the initial conditions. For the above initial conditions 
this is a decreasing function in time. Physically, this means that, unlike the case of an off-resonant oscillation, it takes 
an infinite time for the energy that is transferred from the scissors mode in the xz-plane to the scissors mode in the 
xy-plane to get back to the mode in the a;2:-plane. In Fig. 2 we show this resonance behavior. Finally, we can show 
that in the limit A ^ the oscillatory solution given in Eq. (24) reduces to the nonoscillatory one at resonance given 
by Eq. (30). 



V. SUMMARY AND CONCLUSION 

We have explored the role of the mean-field interaction in coupling the three scissors modes of a Bose-Einstein 
condensate. A variational approach with a gaussian trial wave function, that contains a number of variational 
parameters describing the scissors modes, provides a relatively simple way in which we can extract the main features 
of this coupling. To first order in the deviations in the variational parameters from their equilibrium values we 
reproduce the correct frequencies of the scissors modes. To second order we show that it is not possible to have two 
modes that are coupled if the third mode is not excited. Instead, all three modes need to be excited for nonlinear 
dynamics to occur. 

At higher orders we find that it is possible to couple only two modes. In this case we find a resonance behavior if 
251a;j, — Clxz or Qxy — flxz- Up to the nineth order in the deviations, these are the only two cases of resonance that 
we have found. Close to resonance the equations of motion have been solved exactly using an envelope approach. 
The resulting dynamics is similar to a beating between two modes with a beating frequency 2 A ~ \2Clxy — i~ixz\- We 
notice that the observed resonance behavior [11] occurs exactly for the same condition that we have obtained, namely 
2flxy — ^xz- Also the fact that we did not find any other resonance condition up to the ninth order, indicates that 
the coupling terms that lead to this resonance are also responsible for the experimental resonance in the triaxial case. 

Furthermore, it is important to note here that, quite generally, the down-conversion process from one excitation 
quantum into two excitation quanta with half the energy each, i.e., the so-called Beliaev damping [17], vanishes for 
the scissors modes. This is so because of the negative parity of the scissors modes. In terms of fluctuations of the 
anihilation operator ip{r,t) given by (^(r,t) — '0(r,t) — ■(/'(r, i), the Beliaev damping process is accounted for by an 
interaction term proportional to J dr ■0(r, i)<p^(r, t)(^^(r, t)(^(r, i). Here ip{r,t) =< '(/'(r,t) > is again the condensate 
wave function. Using our variational gaussian wave function given by Eq. (2), we clearly see that with only two scissors 
modes present this integral vanishes, since the integrant is an odd function. We believe that this result is independent 
of our trial wave function and also true within an exact approach [18]. In our case the seventh-order coupling terms 
in Eqs. (18) and (19) correspond to a quadratic coUisional damping process, where three excitation quanta decay into 
four excitation quanta with half the frequency and one excitation quantum with the same frequency. This is shown 
schematically in Fig. 3. This makes sense physically since this is the lowest-order nonvanishing process if Beliaev 
damping is forbidden and we are forced to apply the interaction term proportional to J dr(p'^ {r, t)(^^(r, t)ifi{r, t)(p{r, t) 
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twice to accomplish the down conversion. Another point worth mentioning here is that in the experiments the 
resonance is studied in most detail for an axially symmetric trap with an anisotropy ratio close to A = vZ, when 
the xy-scissors mode is degenerate with the quadrupole mode, which we have excluded in the present calculation. 
In that case we may argue that one scissors-mode quantum can decay into either another scissors-mode quantum 
and a quadrupole-modc quantum or into two quadrupole-mode quanta. These arc in principle also Beliaev damping 
processes, which are however again forbidden for similar reasons as before. Nevertheless, a good understanding of 
these experiments requires an analysis that includes not only the xy and the xz scissors modes but also the quadrupole 
mode, because of the degeneracy that occures in this case. However, our theory should be directly applicable to the 
experiments with a triaxial trap. We hope that in the future more detailed experiments of this kind will also be 
performed, to make a direct comparison between our theory and experiment possible. 

ACKNOWLEDGEMENTS 

We would like to thank C. J. Foot and P. Drummond for useful discussions. This work is supported by the Stichting 
voor Fundamcntecl Onderzoek der Materie (FOM), which is financially supported by the Nederlandse Organisatie 
voor Wetenschappelijk Onderzoek (NWO). 



[1] D. S. Jin, J. R. Ensher, M. R. Mathews, C. E. Wieman, and E. A. Cornell, Phys. Rev. Lett. 77, 420 (1996). 

[2] M.-O. Mewes, M. R. Andrews, N. J. van Druten, D. S. Durfee, C. G. Townsend, and W. Ketterle, Phys. Rev. Lett. 77, 
988 (1996). 

[3] S. Stringari, Phys. Rev. Lett. 77, 2360 (1996). 

[4] K. G. Singh and D. S. Rokhsar, Phys. Rev. Lett. 77, 1667 (1996). 

[5] M. Edwards, P. A. Ruprecht, K. Burnett, R. J. Dodd, and C. W. Clarck, Phys. Rev. Lett. 77, 1671 (1996). 

[6] Y. Castin and R. Dum, Phys. Rev. Lett. 77, 5315 (1996). 

[7] V. M. Perez-Carcia, H. Michinel, J. L Cirac, M. Lewenstein, and P. Zoller, Phys. Rev. Lett. 77, 5320 (1996). 

[8] D. Guery-Odelin and S. Stringari, Phys. Rev. Lett. 83, 4452 (1999). 

[9] M. O. da C. Pires and E. J. V. de Passes, J. Phys. B: At. Mol. Opt. Phys. 33, 3929 (2000). 
[10] O. M. Marago, S. A. Hopkins, J. Arlt, E. Hodby, G. Hechenblaikner, and C. J. Foot, Phys. Rev. Lett. 84, 2056 (2000). 
[11] E. Hodby, O. M. Marago, G. Hechenblaikner, and C. J. Foot, Phys. Rev. Lett. 86, 2196 (2001). 
[12] O. M. Marago, G. Hechenblaikner, E. Hodby, and C. J. Foot, cond-niat/0101213. 
[13] B. Jackson and E. Zaremba, cond-niat/0105465. 
[14] C. J. Foot, private communication. 

[15] G. Baym and C. J. Pethick, Phys. Rev. Lett. 76, 6 (1996). 

[16] F. Dalfovo, S. Giorgini, L. P. Pitaevskii, and S. Stringari, Rev. Mod. Phys. 71, 463 (1999). 
[17] S. T. Behaev, Sov. Phys. JETP 34, 323 (1958). 
[18] See also G. Hechenblaikner, E. Hodby, S. A. Hopkins, O. M. Marago, and C. J. Foot, cond-mat/0105394. 



10 



FIGURE CAPTIONS 



FIG. 1. The real part of Scxz showing two kinds of osciUation. The one with the larger frequency corresponds to the 
unperturbed scissors mode oscillation with frequency ^xz- The slower oscillation is due to the mean- field coupling between the 
scissors mode in the ajz-plane and the xy-p\&ne. The frequency of this oscillation is 2A = 2\2Q.xy — ^xz\. The Bose-Einstein 
condensate parameters used to make this plot are ui^ ~ ij^y ~ ljJz/V^ = 128 Hz. As in the case of Ref. [10], A = 2.54, and 
A'^ = 10"' of *^Rb atoms. The initial conditions are So = 20 x 7r/180 and 0o = 0.03 x tt/ISO. The latter was taken small to 
show that already such a small perturbation in the a;z-scissors mode is sufficient to initiate a substantial coupling between this 
mode and the a;y-mode. 



FIG. 2. The envelope |/(i)| that determines the energy transferred between the two coupled scissors modes. This set of 
curves shows that when approaching the resonance condition A = |2ila;i/ — 0.xz\ ~ 0, the frequency of the curve decreases until 
it becomes equal to zero on resonance. The inset shows this last behaviour for much larger times. For axially symmetric traps 
the resonance takes place at A = \/7. The parameters used for these plots are the same as those mentioned in Fig. 1. 



FIG. 3. A schematic figure representing three down-conversion processes, a) Beliaev damping where one excitation decays 
into two excitations with half the frequency. It is shown in the text that for two scissors modes this process vanishes, b) 
quadratic coUisional damping in which three excitations decay into one excitation of the same frequency and four excitations 
with half the frequency. According to the present work this is the first nonvanishing down- conversion process, c) Beliaev 
damping due to down conversion from one scissors mode into another scissors mode and a degenerate quadrupole mode, or into 
two quadrupole modes. This is the situation that corresponds to the most detailed experiments of Ref. [11], but the matrix 
elements of these proces again vanish. 
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FIGURES 



APPENDIX A: COUPLING MATRICES AND COEFFICIENTS 



The matrices M'^'''=^*hing ^^^^ ]y[scissors ^^^ ^^^^^ ^^ 



I ZL^ 



M 



breathing 



OJ-c ^ ^y ^ ^ z '■ 



25x5 
2 12 2 

7 5 u)^ 5 CJ^ 5 

I — I — W~ 

2 5 TT 5 {jjy 5 



7 5 oJa: 5 ^y 5 

1 1 1~ 

2 5 TT 5 OJj. 5 



7^'^x 


5 t^^ 


3 


257r 


5aj„ 5 


f 


25775 

7« 


aj„ 5 




.i 1 


s 


8 



25 7r5 



^I* '^ ^ z ^ 



7^cj^^^u;j/^ 
I — I E~ 

2 5 TT 5 cj^ 5 

2 22 

7 5 oJx ^ 

"~I — I 1 F 

2 5 TT 5 ijjy 5 a;^ 5 
2 22 2 

-87 5 CJa; 5 CJy 5 

— I — I IE — 

2 5 TT 5 o;, 5 






~1 1 1 

7 5 ujy 5 UJ^ 5 



2 5 TT 5 UJ^ 5 

'~z — n — r 

7 5 uJy 5 cJg 5 



4 1 18 
2'^ 7r"S"cJa; "5" 
~I 1 IT / 

7 5 oJy 5 OJ^ 5 / 



(Al) 



and 



/ 2 12 

r '■f'^^ uJx~^ ^z 
~^ 1 — I S 



M 



scissors 



75 U;^ 5 U)y^ 

— I — I — W 

2 5 TT 5 (jj 5 



1 — I — I — W 

2 5 TT 5 UJy 5 UJz 5 



-UJLUx 



2 12 2~ 

27 5 LJy 5 uJs 5 



■^ 5 o;^ 5 uj-c +ct;2 

3 2 n — 

27^ ojy ^ u;^ "5" 



j'^'^X 






1 18 2 

"2- 



2 12 12 / 

27 5 u)y 5 OJg 5 / 



(A2) 



The coupUng cocfhcients a, j3 and 77, appearing first in Eqs. (18) and (19), are given by 
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and 
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respectively. 
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Fig. 1, U. Al Khawaja 
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Fig. 2, U. Al Khawaja 
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